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ABSTRACT 

Feedback from active galactic nuclei (AGN) has often been invoked both in simulations and 
in interpreting observations for regulating star formation and quenching cooling flows in mas¬ 
sive galaxies. AGN activity can, however, also over-pressurise the dense star-forming regions 
of galaxies and thus enhance star formation, leading to a positive feedback effect. To under¬ 
stand this pressurisation better, we investigate the effect of an ambient external pressure on 
gas fragmentation and triggering of starburst activity by means of hydrodynamical simula¬ 
tions. We find that moderate levels of over-pressurisation of the galaxy boost the global star 
formation rate of the galaxy by an order of magnitude, turn stable discs unstable, and lead to 
significant fragmentation of the gas content of the galaxy, similar to what is observed in high 
redshift galaxies. 
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1 INTRODUCTION 

Supermassive black holes are found at the centers of most, if not all, 
massive galaxies (e.g., Magorrian et al. 1998; Hu 2008; Kormendy 
et al. 201 1). Throughout cosmic history, they are thought to play an 
important role in regulating the baryonic mass content of massive 
galaxies through feedback from Active Galactic Nuclei (AGN) by 
releasing a fraction of the rest-mass accreted energy back into the 
galactic gas and altering the star formation rate (SFR) in the galaxy. 

AGN can exert either negative or positive feedback on their 
surroundings. The former describes cases where the AGN inhibits 
star formation by heating and dispersing the gas in the galaxy, while 
the latter describes the possibility that an AGN may trigger star for¬ 
mation. Negative AGN feedback can operate in quasar-mode from 
radiation at high accretion rates, or radio-mode from AGN jets at 
predominantly low accretion rates (Churazov et al. 2005; Russell 
et al. 2013). It is still unclear how efficiently AGN feedback de¬ 
livers energy (through heating, e.g., Silk & Rees 1998) and mo¬ 
mentum (through physical pushing, King 2003) to the galaxy’s gas 
and what mode of feedback dominates. Both semi-analytical (e.g., 
Croton et al. 2006; Bower et al. 2006) and hydrodynamical cos¬ 
mological simulations (e.g., Di Matteo et al. 2005, 2008; Sijacki 
et al. 2007; Booth & Schaye 2009; Dubois et al. 2010) have shown 
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that negative AGN feedback is an important ingredient in the for¬ 
mation and evolution of massive galaxies, in particular in shaping 
the observed high-end tail of the galaxy mass function, and the low 
SFRs in massive galaxies. Moreover, observations show that cool¬ 
ing flows in the hot circumgalactic and intracluster media can be 
suppressed by the energy transferred by AGN jets (Birzan et al. 
2004; Dunn et al. 2005), again negatively impacting star formation. 

Although AGN feedback has been extensively studied in ob¬ 
servations and through cosmological simulations, the impact on the 
host galaxy and the precise mechanism of the communication of the 
AGN with the galaxy’s interstellar medium (ISM) is far from be¬ 
ing understood. It is not clear why jet feedback, which is thought 
to heat cold gas, should have a similar effect on the multi-phase 
ISM. It has been argued that a jet that propagates through an in¬ 
homogeneous ISM may also trigger or enhance star formation in 
a galaxy (i.e., positive feedback). Begelman & Cioffi (1989) and 
Rees (1989) proposed that the radio jet activity triggers star for¬ 
mation and might serve as an explanation for the alignment of ra¬ 
dio and optical structures in high redshift radio galaxies. Radio jet- 
induced star formation has also been considered as a source pow¬ 
ering luminous starbursts (Silk 2005). Ishibashi & Fabian (2012) 
provide a theoretical framework linking AGN feedback triggering 
of star formation in the host galaxy to the oversized evolution of 
massive galaxies over cosmic time. Furthermore, negative and pos¬ 
itive feedback are not necessarily contradictory (Silk 2013; Zubo- 
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vas et al. 2013; Zinn et al. 2013; Cresci et al. 2015): AGN activity 
may both quench and induce star formation in different parts of the 
host galaxy and on different time-scales. 

Observationally, this positive feedback scenario is directly 
supported by only a few local (Croft et al. 2006; Inskip et al. 2008; 
Salome et al. 2015) and high redshift (Dey et al. 1997; Bicknell 
et al. 2000; Rauch et al. 2013) observations. There are, however, 
also indirect links between jets and star formation which suggest 
possible positive feedback from AGN (Klamer et al. 2004; Mc¬ 
Carthy et al. 1991; McCarthy 1993; Balmaverde et al. 2008; Podi- 
gachoski et al. 2015; Swinbank et al. 2015). 

More recently, high resolution hydrodynamical simulations of 
a jet including a multi-phase ISM have become feasible (Sutherland 
& Bicknell 2007; Antonuccio-Delogu & Silk 2008, 2010; Wagner 
& Bicknell 2011; Gaibler et al. 2012). These studies have shown 
that a clumpy interstellar structure results in a different interac¬ 
tion between the jet and the gas than was assumed from simula¬ 
tions with a homogeneous ISM. It can be generally noted that an 
inhomogeneous ISM affects not only the jet evolution, but also 
the morphology of the host galaxy itself. Simulations by Tortora 
et al. (2009) extended the studies and generalised the simulations 
of Antonuccio-Delogu & Silk (2008) studying the interaction of 
a powerful jet in 2D, two-phase ISM. Tortora et al. (2009) have 
shown that star formation can initially be slightly increased (10- 
20 per cent) followed by a much stronger quenching (more than 
50 per cent) within a time-scale of a few million years. They ar¬ 
gue that the rapid decrease of the SFR after its initial enhance¬ 
ment is a consequence of both the high temperatures as well as 
the reduced cloud mass once the jet cocoon has propagated within 
the medium. Kelvin-Helmholtz instabilities reduce the mass of the 
clouds and, assuming a Schmidt-Kennicutt law, thereby reduce the 
SFR. It should, however, be noted that the 2D approach results in 
a very different temperature and pressure evolution as compared 
to a 3D simulation. Wagner & Bicknell (2011) studied the inter¬ 
action of a relativistic jet interacting with a two-phase ISM at the 
galaxy’s center with a resolution of one kiloparsec. They found that 
the transfer of energy and momentum to the ISM may inhibit star 
formation through the dispersal of gas, but their simulations do not 
contain a star-formation model. It could be argued, though, that be¬ 
cause of the short (< 1 Myr) simulation timescale, the impact of 
cooling is very weak and they therefore might underestimate the 
SFR. 

Gaibler et al. (2012) simulated a powerful AGN jet within a 
massive gaseous, clumpy disc (however they neglected gravity). 
They showed that the jet activity causes a significant change of the 
SFR by enhancing the formation of stars, with inside-out propaga¬ 
tion in the galaxy (see Dugan et al. 2014). Their simulations show 
the formation of a blast wave in the central region of the disc possi¬ 
bly dominating the early evolution between the jet and the galaxy. 
The blast wave results in the formation of a cavity in the disc center 
pushing the gas outwards and compressing the gas within the disc 
at the cavity boundary, generating rings of compressed gas within 
the disc. At later times, the ISM of the disc is pressurised by the 
bow shock enclosing the entire disc. This pressurisation is both due 
to the ram pressure from the backflow but also the thermal pressure 
of the cocoon. It is expected that the thermal pressure dominates 
somewhat as the turbulence is measured to be in the subsonic or 
transonic regime. 

Although the physical understanding of star formation is still 
limited and debated (Padoan & Nordlund 2011), it can be assumed 
that a pressurised disc can trigger gravitational instabilities, com¬ 
press the galaxy’s clouds, and push the densities within the disc 


Table 1. Galaxy parameters: scale radius (r s ), gas fraction (f g ), total stellar 
mass (M*), and total gas mass in the disc (M gas ) 


Identifier 

r s 

fa 

M* 

Mg as 


[kpc] 

[%] 

[10 9 M 0 ] 

[10 9 M 0 ] 

gasLow 

3.4 

10 

8.1 

0.9 

gasHigh 

3.4 

50 

4.6 

4.4 


above the critical density for star formation, thus resulting in an 
increased SFR. This picture is further supported by a few obser¬ 
vations of well-resolved star-forming molecular clouds (Keto et al. 
2005; Rosolowsky & Blitz 2005) and by detailed simulation of the 
ISM (e.g. Slyz et al. 2005; Zubovas et al. 2014). 

Motivated by the pressurisation of the disc found by Gaibler 
et al. (2012) in their hydrodynamical simulations of AGN jet feed¬ 
back, we have investigated the effects of this extra pressure on 
a galaxy disc, by running hydrodynamical simulations with self¬ 
gravity, without AGN jets, but with simple prescriptions for exter¬ 
nal pressure such as that which may be caused by the jet cocoon. 
In a first study (Bieri et al. 2015), we simulated disc galaxies of 
one-tenth the total mass of the Milky Way, varying their initial gas 
fraction. We found that with a given level of external pressure, the 
disk fragments into numerous clumps, causing enhanced star for¬ 
mation. In the present article, we study the effects of external pres¬ 
sure in more detail, by considering different geometries and levels 
of external pressure, as well as studying the effects of supernova 
feedback and mass resolution. 

In Section 2, we describe our suite of hydrodynamical simu¬ 
lations. Our results are presented in Section 3 and summarised in 
Section 4. 

2 SIMULATION SET-UP 
2.1 Basic simulation scheme 

Our simulations begin with a galaxy made of a disc of gas and 
stars, a stellar bulge and a dark matter (DM) halo. We allow this 
galaxy to relax to an equilibrium configuration (with a reasonable 
disc thickness) over the rotation time of the disc at its half-mass 
radius. This first phase is performed without gas cooling, or star 
formation or feedback, in order to evacuate spurious waves emitted 
from the imperfect equilibrium of the initial conditions. After this 
first relaxation phase, we turn on the external pressure, gas cool¬ 
ing, star formation, and also feedback from supernovae (SNe), as 
described below. 

The initial condition method introduced by Springel & Hern- 
quist (2005) is used to generate the DM particles with an 
NFW (Navarro et al. 1997) density profile and a concentration pa¬ 
rameter of c = 10. The virial velocity of the DM particles is set 
to be U 200 = 70kms _1 , which corresponds to a virial radius of 
R 200 ~ 96kpc and a virial mass of M 200 ~ 1.1 x 10 11 M©. 
A Hubble constant of Ho — 73kms _1 Mpc -1 is assumed. The 
star particles as well as the gas are distributed in an exponential 
disc with a scale length of 3.44 kpc and scale height 0.2 kpc, and 
a spherical, non rotating bulge with a Hernquist profile (Hernquist 
1990) of scale radius 0.2 kpc. We use 10 6 DM particles with a mass 
resolution of 1.23 x 10 5 M 0 to sample the dark matter halo, and 
5.625 x 10 5 star particles sampling the disc of which 6.25 x 10 4 
star particles are used to sample the bulge. The stellar mass resolu¬ 
tion is 1.57 x 10 4 M 0 for the 10% gas fraction simulation (here¬ 
after, gasLow ) whereas for the 50% gas fraction simulation (here- 





after, gasHigh) the mass resolution is 8.73 x 10 3 M©. The relevant 
galaxy parameters are shown in Table 1 . 

The simulations are run with the Ramses adaptive mesh 
refinement code (Teyssier 2002). Particles motions are evolved 
through the gravitational force with an adaptive particle mesh 
solver using a cloud-in-cell interpolation, together with the mass 
contribution of the gas component. The evolution of the gas is 
followed with a second-order unsplit Godunov scheme. We use 
the HLLC Riemann solver (Toro et al. 1994) with MinMod total 
variation diminishing scheme to reconstruct the interpolated vari¬ 
ables from their cell-centred values. The box size is 655 kpc with 
a coarse level of 7, and a maximum level of 14 corresponding to 
a Ax = 40 pc minimum cell size for most of the simulations. For 
convergence studies, we perform a higher resolution run with a spa¬ 
tial resolution of Ax — 10 pc (maximum level of refinement 16). 
The refinement is triggered with a quasi-Lagrangian criterion: if the 
gas mass within a cell is larger than 8 x 10 7 M© or if more than 8 
DM particles are within the cell a new refinement level is triggered. 

The circumgalactic medium is modelled with a constant hy¬ 
drogen number density of ticgm = 10 -3 Hem -3 . The pressure 
and temperature profiles outside the disc are calculated assuming 
spherical hydrostatic equilibrium. For the relaxation phase, the sim¬ 
ulations are run for one rotation period of the half-baryonic mass 
radius (5 kpc) of the galaxy, i.e. ^0.5 Gyr. 

The simulations include sub-grid models for cooling, star for¬ 
mation, as well as SN feedback in a subset of runs. The cooling 
mechanism is that described by Sutherland & Dopita (1993), which 
accounts for H, He, and metal contributions to gas cooling (assum¬ 
ing a solar chemical composition of the various metal elements, but 
with a varying metallicity of the gas). The disc is initialised with 
a uniform solar metallicity. No metals are initially placed outside 
the disc. The boundary of the disc is defined using a geometrical 
criteria with cylindrical symmetry using the initial disc radius and 
disc height. Metals are passively advected with gas in the simula¬ 
tion and are modified by individual SNe events with a yield of 0.1, 
which also distribute the metals throughout and outside the galaxy. 
In dense and cold regions, gas is turned into star particles following 
a Schmidt law: 

p*=e*— - it 7lgas > no , (1) 

tfi 

where p* is the star formation rate density, p gas is the gas mass 
density, e* = 0.01 is the star formation efficiency, ts is the local 
gas free-fall time, and n gas and no = 14 H cm -3 for Ax = 40 pc 
(no = 224 H cm -3 for Ax = 10 pc) are the local H number den¬ 
sity and H number density threshold for star formation respectively. 
The Schmidt law is used to draw a probability to form a star with 
a stellar mass of rn* = poAx 3 ~3x 10 4 M© for the low resolu¬ 
tion ( lowRes ) runs and a stellar mass of ?n* ~ 7 x 10 3 M 0 for the 
high resolution ( highRes ) runs (Rasera & Teyssier 2006). The gas 
temperature in high gas density regions (n gas > no) is artificially 
enhanced by a poly tropic equation of state T — To(n gas /no) K-1 , 
where k = 2 is the poly tropic index, and T 0 » 270 K for the 
low and high resolution runs. It is chosen in order to get a constant 
Jeans length resolved with at least 4 cells. This artificial poly tropic 
equation of state is used to prevent the catastrophic and artificial 
collapse of the self-gravitating gas (Truelove et al. 1997). 

We account for the mass and energy release from type II SNe. 
The energy injection, which is purely thermal, corresponds to 

£sN=r?SN^10 5O erg, (2) 

M© 

where psn = 0.2 is the mass fraction of stars going SNe and ra* 



Figure 1 . Mean pressure versus radius at different times (see legend) for the 
pa3 run of the gas High fb set. The dashed line shows the pressure profile 
before the onset of external pressure. The pressures are averaged within 
spherical shells. 


is the mass of the star particle. We also return an amount r]sNm* 
back into the gas for each SN explosion which occurs 10 Myr af¬ 
ter the birth of the star particle. To avoid excessive cooling of the 
gas due to our inability to capture the different phases of the SN 
bubble expansion, we use the delayed cooling approach introduced 
in Teyssier et al. (2013) (in the same spirit as Stinson et al. 2006). 
The energy of the SN explosion is injected into a passive scalar 
variable and blocks the cooling of the gas if the corresponding ve¬ 
locity dispersion is larger than cr t hres = 60kms -1 . The energy 
within that passive scalar decays with a characteristic time-scale 
of tdiss — 2 Myr (tdiss = 0.5 Myr) for Ax m 40 pc resolution 
(Ax = 10 pc respectively), long enough to block the cooling over a 
few cell sound crossing times (see Appendix of Dubois et al. 2015). 


2.2 Application of external pressure 


After adiabatic relaxation (no gas cooling), the origin of time is re¬ 
set to 0 and the base simulations are run further in time with the 
subgrid modeling of gas cooling, star formation and feedback, and 
with an enhanced and uniform pressure outside the disc {pressure 
simulations) for another « 420 Myr. The pressure enhancement 
is applied at an instant starting at t = 0. This instant pressure in¬ 
crease is justified since the bow shock observed in the simulations 
of Gaibler et al. (2012) manages to pressurise the entire gaseous 
disc within a time-frame of only a few Myr. The pressure is en¬ 
hanced in two different configurations: either outside the sphere 
of radius n = 12 kpc (hereafter, pjspher ), or else where the gas 
number density is lower than 0.014Hem -3 (i.e. right outside the 
disc component, hereafter, p_densox equivalently gasHigh jd). In 
the pspher simulation, the bow shock pressurising the disc is as¬ 
sumed to be quasi-isotropic. The effect of isotropy of external pres¬ 
sure is compared with a simulation of a non-isotropic bipolar pres¬ 
sure increase in Appendix A. 

In the case of spherical geometry {pjspher), this pressure en¬ 
hancement is calculated by 


P(r, t) = < 


P(r, 0 ) 

pa 

\r 2 - V! 

pa Pm ax 


r <r i , 

7*1 < r < 7*2 , 

r > r 2 , 


(3) 


where time t = 0 is just before the pressure enhancement, pa is the 
pressure amplification , f{x) — 6x 5 — 15x 4 + 10x 3 is an increas¬ 
ing function of x starting very gradually at x = 0 and smoothly 
reaching a plateau of unity at x = 1, and finally P ma x is the maxi¬ 
mum pressure in the disk at t = 0 (reached in the central few cells), 








with Pmax — 9.8 x 10“ 13 Pa for the gas Low simulation set and 
P max — 4.7 x 10“ 12 Pa for the gasHigh and gasHigh-d simulation 
sets. Here, we adopt 7*2 = r± + 3 kpc. This gradual pressure am¬ 
plification with radius is used to smoothly connect the two pressure 
regimes. For convenience, we will call paX a simulation run where 
the pressure amplification is pa = X. This pressure bath is main¬ 
tained throughout the simulation evolution and is a minimum to 
the pressure evolved in that region. If the pressure within that bath 
becomes larger than paP ma x (due to SNe winds for instance) we 
take the new value of pressure provided by the Riemann solver. For 
this pspher case (but also for the pjdens case), the simulation of no 
pressure amplification corresponds to pa = P(n, 0)/P ma x — 0.1, 
and we will hereafter denote it as nP (for no pressure enhancement). 

For the case of external pressure in disc geometry {pjdens), 
we increase the pressure, only at time t — 0, at a value of pa Pmax 
wherever the gas density is below 0.014 H cm -3 . This gas density 
corresponds to a height of 1.1 kpc above and below the plane along 
the minor axis of the disk (R = 0). 

In the simulations of Gaibler et al. (2012), the bow shock 
that pressurises the disc reaches a maximum pressure of P ~ 
8 x 10“ 11 Pa. This justifies our chosen pressure enhancement 
where the maximum pressure increase for the gasHigh (pa 10) and 
gasLow (pa7) simulation corresponds toP~ 9.8 x 10 -12 Pa and 
P~3.2xl0 -11 Pa, respectively. 

The pressure profiles for one of the pspher simulations (run 
pa3) before and after the pressure enhancement as well as its evo¬ 
lution over the simulation time are shown in Fig. 1. We can see that 
that at 2.4 Myr, right after the restart of the simulation, the pressure 
smoothly rises from ~ 10 -13 Pa at the centre up to 10 -11 Pa at a 
distance r = 13.5 kpc. At later times, this pressure enhancement 
propagates within the central region of the halo and connects to the 
galaxy. 

The relevant physical parameters for the pressure simulations 
are summarised in Table 2. 


3 RESULTS 

In this section, we present our simulation results considering differ¬ 
ent isolated disc simulations with various pressure boosts outside 
the disc. We analyse our simulations regarding disc fragmentation, 
star formation, clump properties, and the galaxy’s mass budget. We 
then compare our simulations with a simple theoretical implemen¬ 
tation regarding the growth of the star formation rate and show that 
it scales approximately as the square root of the external pressure. 
Finally, we calculate the Kennicutt-Schmidt (KS) relation and find 
that our toy model for AGN-induced over-pressurisation leads to 
the galaxies lying higher in the starburst region of the KS relation. 

The effects of external pressure turn out to be similar whether 
or not SN feedback is included in the simulations. We will there¬ 
fore only present, in this section, the results of the stellar feedback 
simulations. A comparison between the non-feedback and feedback 
simulations is provided in Appendix B. 

3.1 Qualitative differences 

Figs. 2 and 3 show maps of the gas density for selected runs 
at different times, for the gasLow Jh and gasHigh Jh as well as 
gasHigh_djh runs respectively, and for two cases without external 
pressure boost nP and with extra pressure pa3. Comparing the nP 
runs, we observe that the gas is clumpier in the gasHigh Jh simula¬ 
tion than in the gasLow Jh run. We will show in Sect. 3.2 that this 


Table 2. Physical parameters of runs: gas fraction (0.1 for gasLow and 0.5 
for gasHigh ), pressure amplification (pa), run with no feedback (nf), run 
with feedback (fb), spatial resolution ( Ax ), and pressure geometry 


Identifier 

gas 

pa 

nf/fb 

Ax 

geometry 


fraction 



[pc] 


paOl = nP 


0.1 

/// 

40 


pa04 


0.4 

/// 

40 


pa08 


0.8 

/// 

40 


pal.2 


1.2 

/// 

40 


pal.5 


1.5 

/// 

40 


pa3 

gasLow 

3 

/// 

40 

pjspher 

pa5 


5 

/// 

40 


pa7 


7 

/// 

40 


paOl = nP 


0.1 

/// 

40 


pa01_hR = nPJiR 

0.1 

x// 

10 


pa02 


0.2 

/// 

40 


pa04 


0.4 

/// 

40 


pa08 


0.8 

/// 

40 


pal.2 


1.2 

/// 

40 


pal.5 


1.5 

/// 

40 


pa2 

gasHigh 

2 

/// 

40 

pjspher 

pa3 


3 

/// 

40 


pa3_hR 


3 

x// 

10 


pa5 


5 

/// 

40 


pa7 


7 

/// 

40 


pa8 


7 

/// 

40 


palO 


10 

/// 

40 


pa01_d = nP 


0.1 

/// 

40 


pa02_d 


0.2 

/// 

40 


pa03_d 


0.3 

/// 

40 


pa04_d 


0.4 

/// 

40 


pa08_d 


0.8 

/// 

40 


pal.2_d 


1.2 

/// 

40 


pal.5_d 


1.5 

/// 

40 


pa2_d 

gasHigh 

2 

/// 

40 

pjdens 

pa3_d 


3 

/// 

40 


pa5_d 


5 

/// 

40 


pa7_d 


7 

/// 

40 


pal0 d 


10 

/// 

40 



is a simple consequence of the Toomre instability. The increased 
pressure leads to accelerated clump formation for the gasLow Jh, 
gasHigh Jh, gasHigh _djh simulations, and a clumpier ISM in all 
cases. Generally less gas between clumps in the enhanced pressure 
runs, in all the gasLow Jh, gasHigh Jh, and gasHigh-dJh cases can 
be seen. 

In Fig. 3, one can compare the two different ways to in¬ 
crease the pressure (pspher, pjdens). The morphological structure 
of the two simulations gasHigh Jh and gasHigh jdjh is slightly dif¬ 
ferent. Fewer clumps are seen in the gasHighJh run than in the 
gasHigh-dJh run. It seems, however, that the clumps are only miss¬ 
ing in the outskirts of the gasHighJh galaxy, whereas a similar 
amount of clumps can be detected in the centre. The edge-on views 
indicate that in the gasHighJh simulations, a large amount of mass 
flows out of the galaxy due to the pressure increase, while in the 
gasHigh-dJh simulations the mass outflow seems to be less ex¬ 
tended. We will quantify the mass outflows in the different runs in 
Sect. 3.5. 

3.2 Disc fragmentation 

Since the star formation recipe depends on the local gas density 
(see eq. 1), we expect enhanced star formation when more clumps 









Figure 2. Gas density maps (mass-weighted) of two of the gas Low jb simulations without enhancement of the external pressure nP (top row), and with 
enhancement of the external pressure pa3 (bottom row). The different columns show different times as labelled. Each panel shows both face-on (40x40 kpc, 
upper part) and edge-on (40x20 kpc, lower part) views. One can see that an increased pressure outside the galaxy leads to accelerated clump formation and 
less gas between the clumps. 


are formed (as the gas gets more concentrated), assuming that the 
clumps have sufficient mass. Therefore, if an increased pressure 
leads to increased fragmentation and hence increased clump for¬ 
mation, we expect star formation to be positively enhanced when 
external pressure is applied to the galaxy. We first consider the 
fragmentation by counting the high-density clumps. We detect the 
clumps in the simulation by running the clump finder described 
by Bleuler & Teyssier (2014). This method identifies all peaks and 
their highest saddle points above a given threshold (21 H cm -3 ). A 
clump is recognised as an individual entity when the peak-to-saddle 
ratio is greater than 1.5; otherwise the density peak is merged with 
the neighbor peak with which it shares the highest saddle point. 

The visual impression of increased clump formation when ex¬ 
ternal pressure is applied on the galaxy (Figs. 2 and 3), is confirmed 
when looking at the number of clumps as a function of time. Fig. 4 
shows the number of clumps as a function of time for the gasLowjb 
(left panel), gasHighjb (middle panel), and gasHigh-djb (right 
panel) simulations, respectively. 

In the gasLowjb run, the number of clumps is constantly in¬ 
creasing with time, regardless of the amount of external pressure. 
Clump formation starts earlier in the runs with external pressure. 
However, the number of clumps at a given time is not a monotonic 
function of external pressure: at low external pressure (up to pa3), 
the number of clumps at given time increases with increasing pres¬ 


sure, while the reverse trend occurs for external pressures above 

3 Pmax. 

The general effect that more clumps are formed in the sim¬ 
ulations with external pressure is similar for the gasHighjb and 
gasHigh-djb runs. Similar to the gasLowjb simulation, the num¬ 
ber of clumps increases with increasing pressure up to a certain 
pressure (pa5) and then decreases again for the gasHighjb simu¬ 
lation and stays at the same level for the gasHigh-djb simulation. 
For the lower pressure as well as the non-pressure simulations, the 
number of clumps increases with time. However, for higher pres¬ 
sure simulations, the number of clumps reaches a plateau at late 
times. For the gasHighjb runs, the initial rise in the number of 
clumps is fastest for the pa3, pa5 and pa7 cases, but in the pa3 case 
the number of clumps reaches its plateau at a later time, hence at a 
higher level. 

The time evolution of the number of clumps for the 
gasHigh-djb simulation is roughly similar to the gasHighjb sim¬ 
ulation. At early times, the rise in number of clumps is fastest for 
the high pressure runs. The number of clumps keep rising with 
time for the lower pressure enhancements, while it reaches a max¬ 
imum for the higher pressure enhancements. The time when the 
number of clumps reaches its plateau is also shortest for higher 
external pressures. After 300 Myr, there is no clear trend in num¬ 
ber of clumps versus external pressure for the higher pressure runs. 
The increase of clump number is therefore highly dependent on the 








Figure 3. Gas density map (mass-weighted) for a selection of the gasHighjb and gas High _djb simulations, for no pressure enhancement (top), and for a 
pressure enhancement of pa3 (middle for gasHighjb and bottom for gasHighAJb). The density scale is as in Figs. 2. An increased pressure outside the galaxy 
leads to accelerated clump formation and less gas between the clumps. The morphological structure of the two simulations with two different ways to increase 
the pressure ( ( gasHighjb , and gasHighAJb ) is slighly different, but only in the outskirts of the galaxy. The edge-on view shows a mass outflow for all the 
simulations. 


pressure applied on the galaxy. However, beyond a certain pressure 
enhancement (pa5_d), the number of clumps remains very similar 
for higher pressure runs. After « 300 Myr the number of clumps 
is roughly independent of external pressure for the gasHighjb and 
gasHigh-dJb simulations. 

In the gasLowJb run with no external pressure, only a single 
clump in the entire disc is formed, at late times (350 Myr), when 
the gas has sufficiently collapsed to reach the clump gas density 
threshold. On the contrary, in the gasHighjb run, the number of 
clumps increases up to ~ 35, even without any forcing by the ex¬ 
ternal pressure. 

The difference between the gas-poor and gas-rich galax¬ 


ies, before the external pressure is applied, is that the gaseous 
disc is Toomre-stable against small-scale fragmentation in the 
gas-poor case where the mean Toomre parameter is (Q) = 

(eg k/ ( 7 tG E gas )) m 3.29 > 1, while the gas-rich disc is Toomre- 
unstable with (Q) s= 0.72 < 1 (see Fig. 5). Here, E gas is the sur¬ 
face density, c s is the sound speed, and k is the epicyclic frequency 
(measuring the shear of the rotating disc). Therefore the gasLowJb 
simulations demonstrate that fragmentation of the galactic disc can 
be driven by the forcing of an external pressure, even though the 
disc is initially Toomre-stable. 
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Figure 4. Time evolution of the number of clumps for a selection of simulations with supernova feedback: gasLow Jb (left), gasHighJb (middle), and 
gasHigh jdjb (right) simulations. The lines are smoothed with a Blackman-Harris window with a width of 2^/len(array), where len(array) is the number of 
points. The clumps were extracted with the Bleuler & Teyssier (2014) algorithm, with a density threshold of 21 H cm~ 3 and a peak-to-saddle threshold of 
1.5. The maps show that an increased pressure leads to increased clump formation and the increase in clump number is dependent on the pressure applied onto 
the galaxy. Beyond a certain pressure enhancement, the number of clumps decreases or remains very similar for higher pressure runs. 


3.3 Star formation history 

In Sect. 3.2, we have seen that an increased pressure enhancement 
leads to an increased number of clumps up to a certain time and 
then a typically lower number of clumps thereafter. Since the gas 
density threshold of clump detection is set to be above that for star 
formation, one expects that the star formation history should evolve 
in a similar fashion to the evolution of the number of clumps. 

Fig. 6 shows that the star formation histories of the different 
runs indeed resemble the time evolution of the number of clumps 
previously shown in Fig. 4. In particular, at early times in the runs 
with gasHigh higher pressures lead to higher SFRs. But with high 
pressures, the SFR saturates earlier. In the gasHighJb runs, the 



Figure 5. Local Toomre parameter for the relaxed disc at t = 0 for the 
gasLow and gasHigh simulations. The dashed line shows the mean Toomre 
parameter of the disc for the gasLow simulation (Q) = 3.29 > 1 and 
for the gasHigh simulation (Q) = 0.72 < 1. One therefore expects 
the gasHigh simulation to fragment independently of external pressure en¬ 
hancement whereas the gasLow simulations are not expected to fragment 
into many clumps. The fact that the pressur-enhanced simulations of the 
gasLow galaxies shows a significant increase in the number of clumps com¬ 
pared to no-pressure enhancement demonstrates that external pressure can 
stimulate the fragmentation of a disc even if it is Toomre-stable. 


maximum SFRs in the high pressure runs are lower than in the other 
runs, while in the gasHigh-dJb runs, the maximum level of SFR 
is reached for the three highest pressures, while the SFRs at later 
times (300 Myr) are roughly independent of the external pressure. 

In the gasLow Jb runs, while the nP case leads to star forma¬ 
tion only after a long time delay (330 Myr), the highest pressures, 
although leading to immediate but small levels of star formation, 
are unable to generate substantial star formation from the earliest 
times. The runs with intermediate pressures produce the highest 
SFR at all times. We will show in the next sections that this is due 
to the low mass-outflow of the intermediate pressure simulations 
that allows the clumps to increase in density. Conversely, larger 
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Figure 6. SFR for a selection of the simulations with supernova feedback: gasLow jb (left), gasHighjb (middle), and gasHigh-djb (right) simulations. The 
lines are smoothed as in Figure 4. This figure shows that higher pressures lead to higher SFRs. For the higher pressure enhancement simulations, the SFR 
reaches a maximum or plateau at later times for the gasHighjb, and gasHigh-djb simulations. After a certain pressure increase, the SFR decreases or stays at 
the same level for all the simulations and hence the runs with intermediate pressure generally produce the highest SFR at all times. 


external pressures lead to such strong pressure waves that the gas 
is removed from the galaxy. This prevents the formation of large 
clumps and tends to suppress the star formation. 

The effect of the external pressure on the SFR is even more 
significant when looking at the gasHighjb simulations (middle 
panel of Fig. 6). The SFR of the no-pressure simulation slowly in¬ 
creases after a certain time, whereas the SFR increases faster when 
pressure is applied: it reaches a maximum at a certain rate and more 
or less maintains this rate for the remaining of the simulation. To¬ 
wards the end of the simulation, the SFR of the no-pressure simu¬ 
lation catches up, again similarly to the clump number behaviour. 
The SFR for the gasHigh-djb simulation (right panel of Fig. 6) 
behaves quantitatively similar to the gasHighjb simulation. One 
can see that the SFR in these simulations reaches the maximum or 
plateau at later times than in the corresponding gasHighjb simu¬ 
lations. The rapid rise of the SFR reaches increasingly higher lev¬ 
els of peak SFR with higher external pressure up to pa5_d, while 
palCLd reaches a slightly lower maximum SFR. 

The left panels of Fig. 6 show that the SFRs of the gasLow 
simulations start with a significant time delay, and the maximum 
enhancement of the SFR relative to the nP run is highest (~ 12) 
at the end of the simulation (after 400 Myr). On the other hand, 
the corresponding SFR enhancements for the higher gas fraction 
simulations (middle and right panels of Fig. 6) are lower (~ 3.5 
for gasHigh and ~ 1.5 for gasHigh-d) at the end of the simulation 
(after 400 Myr) than at the beginning (~ 40 for gasHigh and ~ 70 
for gasHighud at ~ 80 Myr) of the simulation. External pressure 
thus first produces a significantly higher SFR in comparison to the 
simulation with no external pressure. But the duration of this large 
SFR enhancement for the gasHigh and gasHigh-d simulations is 
shorter than that of the gasLow simulation. The free fall time of 
the higher density gas is shorter than the free fall time of the low 
density gas which leads to the gas collapsing early on, whereas a 
delay is expected for the lower gas fraction disc. 

3.4 Clump Properties 

An important SFR requires a significant supply of cold gas as well 
as the fragmentation of the disc into clumps that carry a sufficient 
amount of gas to form stars. On the other hand, one can argue with 


the Jeans and Toomre instability arguments if indeed an increased 
pressure outside the galaxy that later increases the pressure inside 
the galaxy leads to higher densities, as well as a possible expulsion 
of disc gas depending on the momentum carried by the pressure 
wave coming into the disc. The competition between higher densi¬ 
ties and mass outflow will influence the amount of gas within the 
clumps. 

For the gas-rich disc simulations, we saw (Fig. 4) that, at the 
very beginning, when the pressure wave comes into the galaxy, the 
number of clumps is highest for the highest pressure. While the 
clumps are more numerous with the highest pressures, it is worth¬ 
while knowing whether their masses are affected by the external 
pressure. 

Fig. 7 shows the modulation of the average clump mass rises 
with external pressure at three times of the simulations. At early 
times (top panel), the clump masses for the gasHigh simulations 
{gasHigh, center panel, and gasHigh_d, right panel) are higher the 
greater the external pressure. At later times (middle and bottom 
panel), the clump masses are roughly independent of the external 
pressure applied, probably because the disc gas has been either al¬ 
ready accreted onto the clumps or expelled out of the galaxy (see 
discussion below), leaving no more diffuse gas available for accre¬ 
tion onto the clumps. The time at which the diffuse gas is either 
consumed onto the clumps or expelled must happen earlier for the 
higher external pressure simulations as the fragmentation happened 
earlier for these simulations. This explains the different times when 
the SFR reaches a plateau, occurring earlier the higher the pressure. 
In the gasHighjb runs at high external pressures, fragmentation is 
not the only cause of SFR (since there is a maximum pressure en¬ 
hancement beyond which the SFR is lower), meaning that the gas 
supply is more critical, and not always available despite the high 
gas fraction. This suggests that the mass flow out of the galaxy also 
plays an important role. And indeed, the mass outflow is very effi¬ 
cient for pa7 and palO after 30 Myr. We will discuss this in detail 
in Sect. 3.5, below. 

In the p-dens simulations, there is a maximum pressure en¬ 
hancement (pa5_d) beyond which the SFR remains at approx¬ 
imately the same level without decreasing. As we will see in 
Sect. 3.5, mass outflows are also absent. The high gas fraction 
leads therefore to higher density enhancement by external pressure, 
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Figure 7. Average clump mass for the gasLow Jb (left), gasHighJb (middle), and gasHigh A Jb (right) simulations. The orange line with round markers 
correspond to the non-feedback simulations whereas the pink line with triangles corresponds to the feedback simulations. The shaded green/blue area shows 
the area of mass containing Per of the density PDF for the no-feedback/feedback simulations. In the bottom and top of each subfigure, the x-axis shows the 
P j Pm ax and P/P 20 kpc values respectively, where P max is the maximum pressure inside the disk and P 2 o k P c is the averaged pressure at 20 kpc. At the 
beginning of the simulation, the clump masses for the gasHigh simulations are higher the greater the pressure for the gasHigh simulations. At later times, the 
clump masses are roughly independent of the external pressure. For the gasLow simulation, the clump mass does not increase with higher external pressure 
but rather decreases or stays at approximately the same level. 
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Figure 8. Density PDF at different times for a selection of the gasLow Jb (left), gasHigh Jb (middle) and gasHigh A Jb (right) simulations. The threshold 
(14Hem -3 ) for star formation is plotted as a grey vertical line. One can see that greater external pressure allows the galaxy to reach higher densities on 
a faster time-scale. The nP simulation slowly catches up with the over-pressure simulations for the gasHigh whereas this is not the case for the gasLow 
simulation. For the gasLow simulation, the high pressure simulation (pa7) never reaches the high densities. As indicated with the grey power law lines one 
can see that over the course of the simulation a slope between -3/2 and -7/5 (or even steeper) develops at high densities. This is in agreement with simulations 
including gravity and turbulence done for instance by Kritsuk et al. (201 1). 





































hence both number of clumps and SFR are highest when the exter¬ 
nal pressure are high. However, there only is a limited amount of 
gas available in the galaxy. One can assume that the limited gas 
supply is insufficient for more star formation, so the SFR remains 
at the same level independently of the pressure enhancement. 

Fig. 7 shows that, for the gasLow simulations ( gas Low , left 
panel), the clump mass does not increase with greater external pres¬ 
sure, but rather decreases at early times (top panel). At later times, 
the clump mass is roughly independent of the pressure up to pa3, 
beyond which the clump mass decreases. This is the same pressure 
enhancement which leads to the highest SFR. As we will see in 
Sect. 3.5, the gasLow galaxies suffer from strong gas outflows that 
reduce the supply of gas available for clump buildup, leading in 
turn to smaller clump masses within the galaxy. 

In Fig. 7, the difference in clump masses for the feedback 
and non-feedback simulations can also be seen for all the simu¬ 
lations. At early times, there is a significant difference between the 
no-feedback and feedback simulations. At later times, the differ¬ 
ence in clump masses becomes more apparent for both the high 
gas fraction and gasLow simulations. One can see that the clump 
masses for the feedback simulations are lower at the end of the 
simulation than for the no-feedback simulations, independent of 
the pressure increase. Because the feedback increases the poros¬ 
ity of the interstellar medium that in turn counteracts the formation 
of clumps (Silk 2001), the observed smaller clump masses for the 
feedback simulations are expected. This difference is more domi¬ 
nant in the gasHigh simulation as can be seen in the left panel of 
Fig. 7. The similar mean clump masses in the feedback and non¬ 
feedback runs at the beginning of the simulation appears to be a 
consequence of the implementation of the SNe in the simulation. 
As discussed above, a SN explosion occurs 10 Myr after the birth 
of the star particle. The first stars form shortly before 50 Myr and 
one would therefore not expect to see a large difference between 
the feedback and no-feedback simulations. At 100 Myr, some stars 
exploded into SNe, but only constitute a small fraction of all stars, 
hence the small difference between the feedback and no-feedback 
simulations at this stage. 

It is interesting to look at the density probability function 
(PDF) at different times of the simulations in order to better under¬ 
stand the observed SFR behaviour. The PDF can be seen in Fig. 8 
for two different times for a selection of the gasLow (left), gasHigh 
(middle), and gasHigh A (right) simulations. Increasing external 
pressure allows one to reach higher gas densities faster, which is 
in agreement with the SFR behaviour we have seen previously. For 
the gasHigh simulations, the nP simulation slowly catches up with 
the over-pressure simulations similar to the SFR behaviour. For the 
gasLow simulation, the no-pressure simulations never attain the 
densities reached by the moderate pressure enhancement simula¬ 
tions. The high pressure simulation (pa7) also never reaches high 
gas discuss in Sect. 3.5, below. It can be seen in Fig. 8 that over 
the course of the simulation a high density power law with a slope 
between -7/4 (or even steeper) and -3/2 develops, especially for the 
high gas fraction simulations. A comparable power law range has 
been found in observations (e.g., Kainulainen et al. 2009, Lombardi 
et al. 2010) and simulations including gravity and turbulence (e.g, 
Kritsuk et al. 2011) where they argue that the origin of the power 
law tail is due to self-similar collapse solutions. 

3.5 The Galaxy’s Mass budget 

The mass flow rate (MFR) as well as the total amount of newly 
formed stars plus dense gas should provide us a better understand¬ 


ing of the star formation history described above. In particular, one 
would like to understand why there seems to be an optimal external 
pressure enhancement for star formation, beyond which the SFR 
ends up at lower values. 

We measure the gas mass flux through a sphere of radius 
16 kpc as 

M gas = §§ p v • r dS = m i v i • ^ , (4) 

iGshell 

where i denotes the index of a cell within a spherical shell of sur¬ 
face S and volume V. Here, we adopt a shell of thickness 4 kpc. 
The MFR is shown in Fig. 9 for the gasLow Jb (left), gasHigh Jb 
(middle), and gasHighjdJb (right) simulations, again only for the 
simulations with SN feedback. The top panels show the mass flow, 
while the bottom ones show the cumulative mass flow in fractions 
of the total gas mass within the 16 kpc sphere before the pressure 
increase shown on the bottom. We will discuss the effects of SN 
feedback in Appendix B. 

In all three sets of simulations, external pressure leads to mass 
inflow at early times. This early mass inflow is large but different 
for the different ways pressure is applied onto the galaxy. In the 
pjspher simulations, the pressure is applied outside the galaxy in 
a low density medium leading the pressure to have a larger pres¬ 
sure gradient than in the pjdens simulations where the pressure is 
applied close to the galaxy and therefore in a higher density envi¬ 
ronment. This larger pressure gradient in the pspher simulations, 
allows the pressure wave to carry more mass and momentum from 
the ambient hot medium in comparison to the pressure wave of the 
P-dens simulations. This explains the larger mass inflow observed 
in the pspher simulations (left and middle panels of Fig. 9) com¬ 
pared to that at the start of the pudens simulation. With its larger 
momentum, the mass inflow of the pjspher pressure wave is fol¬ 
lowed by a short period of mass outflow (for both low and high 
gas fractions). This mass outflow is negligible for the gasHigh-djb 
simulations, since the pressure wave carries very little momentum. 

For the pspher simulation sets, higher external pressures lead 
to stronger maximum inflows at early times and to stronger max¬ 
imum outflows at later times. In addition, in the simulations with 
high external pressures (pa7 and pa 10), the mass outflow that fol¬ 
lows the mass inflow occurs very rapidly (in less than 20 Myr). 
After these two phases of important mass inflow/outflow, the MFR 
oscillates around zero for both the gasLow Jb and gasHigh Jb sim¬ 
ulations with pa < 5. In contrast, in the gasHigh Jtjb simulations, 
the MFR depends little on the external pressure. 

It is also instructive to consider the evolution of the cumula¬ 
tive gas mass flow through the 16 kpc sphere (bottom panels of 
Fig. 9) shown relative to the initial mass within the sphere of that 
radius. The cumulative mass flow remains negative (e.g. mass in¬ 
flow) for all gasLow Jb simulations. The cumulative mass flow for 
the gasHigh Jb with low external pressures (pa0.2 and pa0.8) re¬ 
main negative, while for higher external pressures, they end up 
positive. Whereas for low external pressure (< pa3) the cumula¬ 
tive mass outflow is less than 5% of the initial gas mass, the cu¬ 
mulative mass outflow reaches 80% of the initial gas mass for the 
highest external pressure simulation palO. Finally, all simulations 
in the pjdens geometry lead to cumulative mass inflow at all times, 
with the strongest cumulative inflows occurring for the runs with 
the greatest external pressures. We stress, however, that the pres¬ 
sure and no-pressure pjdens simulations do not differ significantly 
and that, overall, there is little net mass flow. This most likely ex¬ 
plains why the SFR of the p-dens simulations is smoother and less 
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Figure 9. Time evolution of the mass flow rate (top) and cumulative mass flow relative to the initial mass Mo,tot before inputting the pressure (bottom), both 
at 16 kpc from the galaxy center (where the initial mass is within a sphere of that radius), for selected runs from the gas Low Jb (left), gasHighJb (middle) and 
gasHighAJb (right) simulations. Negative (positive) values of the mass outflow rate denote a net mass inflow (outflow). One can see a difference in the mass 
flow rate (MFR) between the different ways pressure is put on the galaxy ( gasHighJb , gasHighAJb). Due to the pressure gradient in the pspher simulations, 
the pressure wave coming into the galaxy carries a lot of momentum that leads to a mass inflow followed by a mass outflow for the higher pressure simulations. 
This mass outflow is negligible for the pjdens simulations due to the pressure wave carrying very little momentum. For the most extreme case, the expelled 
mass reaches 80% of the initial gas mass for the highest external pressure simulation of gasHighJb. 
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Figure 10. Time evolution of the mass in newly formed stars plus dense (n > 0.1 H cm -3 ) relative to the initial gas mass for a selection of the simulations 
with SN feedback: gasLow Jb (left), gasHighJb (middle), and gasHighAJb (right) simulations. The lines are smoothed as in Figure 4. Due to the mass inflow 
shown in Fig. 9, at the beginning of the simulation more mass can end up in the galaxy. This extra mass of gas is more significant for the gasLow galaxy 
for both the no-pressure and low-pressure enhancement simulations. For high pressure enhancement in the gasLow simulation, the incoming pressure wave 
significantly disperses the galactic gas. The evolution of Mtot is similar for the gasHighJb simulation albeit with less mass variation. For the gasHighAJb 
simulations, no significant mass variation within the galaxy due to pressurisation is observed. 


noisy when compared to the pjspher simulations of the same gas 
fraction (gasHigh). 

In order to understand the galaxy’s mass budget better we look 
at the time evolution of the total mass of newly formed stars plus 
dense (ti A 0.1 Hem ) gas, .Mtot — -Mtot,starsN -Mtot,gasD is 


shown relative to the initial galaxy gas mass (total not just dense), 
M to t,gaso (see Table 1) shown in Fig. 10. The initial value of 
Mtot /Mtot,gaso is below unity at t = 0 since the gas density in 
the galaxy is not everywhere above n > 0.1 Hem -3 , especially 
in the outskirts of the disc and for the gasLow galaxy. In the ab- 





























sence of extra external pressure (nP runs), the ratio M to t /M t O t, g aso 
quickly moves significantly beyond unity as the gas cooling allows 
to reach the gas densities above n > 0.1H cm -3 . The gas cooling 
also takes place in the circumgalactic medium that feeds the galaxy 
with some extra gas. This extra mass of gas adds more significantly 
to the low-gas fraction galaxy because of its lower initial gas mass, 
which explains why the increase is more significant in the gasLow 
runs than in the gas High runs. 

For the gasLow simulations, intermediate regimes of forced 
external pressure (from pa0.4 to pa3) also show values of 
M to t/Mtot,gas 0 above 1 with values comparable to the nP run. 
Therefore, the increase in SFR (an order of magnitude above nP) 
is to be attributed to the extra compression of the ISM and explo¬ 
ration of larger gas densities with shorter collapsing time-scales 
(see Fig. 8). In contrast, higher pressurisation values of the ISM 
(pa5 and pa7) lead to strong gas removal due to the incoming pres¬ 
sure wave that manages to significantly disperse the galactic gas. 
Since the gas reservoir is reduced, the SFR is also suppressed com¬ 
pared to more intermediate regimes of pressurisation, but the over¬ 
all SFR is still larger than in the nP case, where gas fragmentation 
is not reached. 

The evolution of Mtot/M to t,gasO in the pspher gasHigh 
galaxy behaves similarly to that for the gasLow galaxy, although 
with lower mass variation. It starts below unity for all pressures, 
and decreases even more for the high pressure increases (pa5 and 
palO), because of the large mass outflows observed for those runs. 
This shows that the large mass outflows associated with high pres¬ 
sures prevent star formation. On the other hand, M to t/Mtot,gas0 
keeps rising for the lower pressure enhancements, showing that be¬ 
cause no large mass outflow is observed, more stars can be formed. 
The M to t/Mtot,gas 0 curve is higher when a small pressure is ap¬ 
plied outside the galaxy compared with the no-pressure simula¬ 
tion. In the case of p_dens over-pressurisation, there is no signif¬ 
icant (< 20 percent relative) mass variation in the galaxy. Thus, 
the early fragmentation due to the increased pressure drives the dif¬ 
ferent SFR levels for the different pressure simulations. 

3.6 The Star Formation Rate 

We expect the SFR to scale as the square root of the external pres¬ 
sure, at least in the gasHigh case, for the following reason. The 
Kennicutt (1998) star formation relation is generically fit by 

s* = (5) 

^dyn 

for theoretical and observational reasons (Elmegreen 1997; Silk 
1997; Genzel et al. 2010). In equation (5), E* is the surface density 
of star formation, E g is the gas surface density, e is a dimensionless 
normalisation constant, and tdyn is the dynamical time (the rotation 
time for a disk galaxy or the free-fall time for a giant molecular 
cloud, both classes of objects fitting the correlation Krumholz et al. 
2012). The slope of the correlation corresponds to the linear theory- 
inspired convolution of gas density and the most rapidly growing 
gravitational instability rate for a cold disk. Much more physics re¬ 
sides in the normalisation, e, which is a measure of the star forma¬ 
tion efficiency, often defined as the fraction of gas turned into stars 
per dynamical time, and in the dispersion. We will not address the 
dispersion here, other than to remark that physics beyond a density 
threshold must be included, as is evident from the low (e.g. cen¬ 
tral molecular zone Kruijssen et al. 2014) and high (both in nearby, 
cf. Leroy et al. 2015; Turner et al. 2015 and distant, cf. Finkelstein 


et al. 2015; Dye et al. 2015) outliers. The usual fit to the normali¬ 
sation is e « 0.02, to star-forming systems out to z ~ 2, although 
there is recent evidence that the efficiency is significantly higher in 
starburst galaxies at higher redshift. The logical generalisation to 
go beyond a density threshold is to include turbulence (Padoan & 
Nordlund 2011; Hennebelle & Chabrier 2011; Hopkins et al. 2013). 
An especially simple implementation is given in Silk (2001); Silk 
& Norman (2009). We will now extend their argumentation of mo¬ 
mentum injection by SNe to an injection of energy by any kind of 
process that affects the whole galaxy (e.g. energy injection given 
by external pressure). 

Let Eiui be the kinetic energy injection given for instance by 
a pressure wave coming into the galaxy. In the same fashion let 
^inj^inj be the momentum injection. We now make the educated 
assumption that the energy injection affects the interstellar clouds. 
These clouds acquire terminal velocity given by 

tcoll ,s\ 

o-g = e I'inj 7- . (6) 

t-dyn 

where t C oii is the cloud collision timescale and is the specific 
momentum injected. We immediately see that the key efficiency 
parameter e is proportional to the gas turbulence velocity a g . There 
are recent indications that star formation efficiency increases in 
highly turbulent environments (Leroy et al. 2015). With momen¬ 
tum balance and Eq. (6), we get 

E* -Einj _ fc Eg CTg 
mini V c tcoll 

where cr g is the gas velocity dispersion, v c the cloud velocity and f c 
the cloud volume filling factor. We however ignore complications 
with f c and define it directly by the cold gas fraction / g . Also, we 
assumed Ui n j oc v c , with v c being the velocity of the cloud. We can 
rewrite the star formation rate per unit volume as 

P* — Cinj fg \/G pg pg (8) 

With €inj — (minj U c CTg)/_E/i n j . 

For the present purpose, we note that the SFR can be rewritten 
as 

P* = — G 1/2 p e pl^ h , (9) 

°g 

where we used p g = p g and assumed that p g oc p t urb, with 
Pturb being the turbulent pressure induced by the injection. 

Hence, we expect that AGN-induced pressure should provide 
a boost of the star formation rate, independently of any possible in¬ 
crease in star formation efficiency, and initially vary as the square 
root of the pressure boost. This effect should be important for the 
gas-rich simulations as gas turbulence plays a lesser role in con¬ 
tributing to a gas-poor disk and the assumption that the gas pres¬ 
sure scales with the gas turbulence breaks down. This is because 
we omitted the cold cloud mass fraction in the simple formulation 
above, and the molecular hydrogen fraction decreases at lower gas 
pressure (Blitz & Rosolowsky 2006). 

Fig. 11 illustrates how the SFR increases in time as a function 
of the external pressure. The mean SFR is calculated from a linear 
fit of the SFR from time t = 0 to time t. For the gasLow , gasHigh , 
and gasHigh-d simulation sets, the SFR increases with increasing 
external pressure until a maximum is reached and then decreases 
again or stays at a similar level as for the gasHigh-d simulations. 
Also, no significant difference in the SFR can be seen between the 
feedback and non-feedback simulations, indicating that the external 
pressure increase is the dominant effect driving the increased SFR. 

In order to compare with the prediction explained above, a 
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Figure 11. Time evolution of the average rate of increase of the star formation rate for all the gasLow Jb (left), gasHigh Jb (middle), and gasHigh A Jb (right) 
simulations. Here, the no-feedback simulations are shown in blue and the feedback simulations are shown in green. The average rate of increase is calculated 
from a linear fit of the SFR from t =0 to three different times. The red points are the corresponding points from the highRes runs. The linear fit has been 
performed on the smoothed SFR data. The SFR data are smoothed as in Figure 4. The dark red-dashed straight lines correspond to visual fits with a slope 
1/2 to guide the eye. One can see that for higher gas fraction simulations, the SFR growth rate increases with external pressure and follows the square root of 
external pressure. At later times, this is only true up until a certain external pressure. For the low gas fraction simulation, the SFR growth rate does not scale 
with the square root of the external pressure. 
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Figure 12. Kennicutt-Schmidt (KS) relation for selected runs (nP, pa3, pa7 with blue, green and red curves, respectively, with lighter colours for the runs 
with 4 times higher spatial resolution) from the gasLow Jb (left), gasHigh Jb (middle) and gasHigh A Jb (right) simulations at 16 different times in color. 
The evolution in time is shown with arrows. The 16 different times are equally spaced between f s tart = 25 Myr and t en d = 400 Myr. The circle shown 
corresponds to t en ^ of the simulation with the corresponding color. Data from Daddi et al. (2010) for the starburst (dashed line) and quiescent (solid line) 
sequences are over-plotted for reference. For all the simulations, independently of the gas fraction, the simulations with external pressure lie closer to the 
starburst sequence than the no-pressure simulations. This shows that our toy model for AGN-induced star-formation might be a possible explanation for 
explaining the increased number of starburst galaxies observed in the distant Universe. 





























dark red dashed curve is plotted to guide the eye in Fig. 11 rep¬ 
resenting a square root fit of the SFR as a function of the exter¬ 
nal pressure applied on the galaxy. One sees that, for the gas-rich 
gasHigh and gasHigh-d simulations, the SFR follows the square 
root of the external pressure very well. At later times, the square 
root law is only fulfilled until the optimal pressure is reached. 
For the gasLow simulations, the SFR does not scale well with the 
square root of the external pressure. But as explained above, this 
was not unexpected. 

The bright red points in Fig. 11 are the corresponding points 
from the highRes run. As can be seen in Appendix C the formation 
of stars in the highRes run shows a delay compared to the lowRes 
run. This is due to the increase in star formation threshold for the 
higher resolution run. Because of this delay, one can also see a de¬ 
layed behaviour of the SFR of the highRes run in Fig. 11. However, 
at later times, the highRes simulation shows a similar behaviour to 
the lowRes simulation. 

3.7 The Kennicutt-Schmidt Relation 

The Kennicutt-Schmidt (KS) law (Kennicutt 1998) law relates the 
SFR per unit area as a power of the surface density of gas. This 
relation holds over several orders of magnitude in both quantities 
(Krumholz, Dekel, & McKee 2012, and references therein), with 
the same normalisation for global galaxies (including high redshift 
ones, Genzel et al. 2010) and giant molecular clouds in the Milky 
Way (Heiderman et al. 2010; Lada, Lombardi, & Alves 2010) and 
in the nearby M51 galaxy (Kennicutt et al. 2007). This demon¬ 
strates the remarkable universality of the SFR. At high redshift, 
starburst galaxies lie above the KS law for normal galaxies (Genzel 
et al. 2010). While the cause of this observed offset is not known, 
one may speculate that this increased SFR may be caused by posi¬ 
tive AGN feedback. 

To investigate this in further detail, we check whether our 
pressurised galactic discs follow the KS relation for normal galax¬ 
ies, or are above it as starbursts are observed to be, or below it. 
We adapt here an equivalent technique to Powell et al. (2013) and 
calculate the half-light radius by assigning a luminosity to each star 
particle dependent on their age and proportional to their mass (Wei- 
dner et al. 2004), 

Rage < 10 Myr) oc M stars (10) 

/ \ -0 ' 7 

Rage >10 Myr) oc M star3 ( J (H) 

We randomly assign an age in the range 0-5 Gyr for stars that are 
specified in the initial conditions and therefore have an age equal to 
zero when the simulation starts. For a given output, Esfr is calcu¬ 
lated within the half-light radius using the SFR averaged over the 
previous 10 Myr. The KS relation is calculated by finding the 3D 
half-light-radius. Within this volume, all the gas above a threshold 
of 0.1 Hem -3 is used to calculate E gas and all the new stars are 
used to calculate Esfr, however the quantities are divided by the 
area i rr| D . 

Fig. 12 shows the KS relation at different times for sev¬ 
eral simulations for the three cases of gasLow Jb, gasHigh Jb, and 
gasHigh-djb with the observed relation from Daddi et al. (2010) 
overplotted. One can see that, over the course of 400 Myr, all runs 
lead to an increase in Esfr, by one to two dex, with much smaller 
variations (less than 0.3 dex) in E gas . In particular, the runs with 
high gas fraction (with or without external pressure) show a de¬ 
crease in E gas . This decrease in gas surface density is related to the 


gas mass outflow at late times (Fig. 9) and to the consumption of 
gas by star formation. For the gasLow simulation, the evolution of 
E gas is tied to the evolution of the total baryonic mass within the 
galaxy shown in Fig. 10. The light blue and green lines in Fig. 12 
show the highRes simulations for the nP and pa3 set, respectively. 
One can see that the trends of the highRes simulation are very simi¬ 
lar to the trends of the lowRes simulations, specifically at late times. 
Especially for the pa3 highRes simulation, one can see that higher 
gas densities are reached due to the higher resolution which allows 
the gas to collapse even further. 

For all the simulations independently of the gas fraction, the 
simulations with external pressure end up being pushed closer to or 
further beyond the starburst sequence than the corresponding simu¬ 
lations without external pressure. This trend is not changed for the 
highRes runs. Runs with external pressures leading to higher SFR 
also have a higher Esfr and therefore end up even closer to the 
starburst sequence. For instance, for the pspher simulations, the 
pa7 run does not extend as far beyond the observed starburst se¬ 
quence as the pa3 run, which indeed reaches higher SFR for both 
the gasLow Jb and gasHigh Jb simulations (Figs. 6 and 11). This 
is not as significant for the gasHigh_dJb simulations as they do 
not experience a decrease of SFR but rather that the SFR stays at a 
certain level after a certain pressure increase (~ pa5). We therefore 
see for this simulation set that the KS relations end in a very similar 
parameter space. 

We conclude that this toy model for AGN-induced over- 
pressurisation plausibly leads to AGN-associated star-forming 
galaxies having enhanced specfic star formation rates, for example 
as suggested by recent observations, cf. (Zinn et al. 2013; Drouart 
et al. 2014). 


4 CONCLUSIONS 

It is a fascinating challenge to understand the extreme star forma¬ 
tion rates observed for some high-redshift galaxies, typically with 
luminous AGN and massive outflows: are these caused by higher 
contents of molecular gas or by a greater efficiency of star forma¬ 
tion relative to this molecular gas content? Is turbulence sufficient 
to explain the high SFR values, or do we need recourse to a more 
exotic pathway that enhances star formation rates even more? The 
latter option is motivated by the increasing evidence for the role 
of AGN in star formation, and in particular their role in a putative 
phase of positive feedback that accompanies or even precedes the 
commonly observed massive, star formation-quenching, outflows 
stimulated by AGN activity. 

Using hydrodynamic al simulations of isolated disc galaxies 
embedded in a hot over-pressurised halo, we have been able to 
study the response of the galaxy SFR to the forcing exerted by this 
external gas pressure onto the disc. The pressure enhancement trig¬ 
gers instabilities leading to more fragmentation when compared to 
the no-pressure simulations (Figs. 2 and 3). The enhanced frag¬ 
mentation leads to the formation of more clumps (Fig. 4) as well as 
larger values of SFR (Fig. 6). This hints at a positive effect of the 
pressurisation of the disc and therefore to positive feedback. 

We observe a difference in the behaviour for the different ways 
in which the pressure is applied. In the simulations where external 
pressure is continuously applied beyond a certain radius (pspher 
simulations), we observe an optimal pressure beyond which the 
number of clumps as well as the SFR is decreased. For the simula¬ 
tions where the pressure is instantaneously applied using a density 



threshold (over-pressure applied closer to the galaxy disc), such an 
optimal pressure is not observed. 

We have seen that the mass outflow plays a role in explaining 
this optimal pressure. In particular, for the gasHigh-fb simulations, 
a significant amount of gas gets expelled out of the galaxy, leaving 
little gas left to form stars and thereby lowering the SFR. The dif¬ 
ference in SFR between the high and low external pressures for the 
gasLowjb simulations is explained by the stagnation of the accu¬ 
mulation of mass in the clumps, which is again related to the large 
amount of gas that is removed by the incoming pressure wave. Our 
simulations have been tested with respect to the resolution and lo¬ 
cal presence or absence of SN explosions: the over-pressurisation 
of the disc still leads to a positive feedback effect (enhanced SFR). 

We found that at given times of the pspher simulations, the 
SFR (and its mean growth rate) vary as the square root of the ap¬ 
plied pressure. We explain this by adapting the Schmidt law for 
the SFR as a function of 3D gas density for the inclusion of ex¬ 
tra pressure caused by the AGN bow shock-driven radio lobe or 
wind, leading to compression times typically an order of magni¬ 
tude shorter than the dynamical time, as argued by Silk & Norman 
(2009). 

Though our setup of the extra pressure exerted by circum- 
galactic gas onto the galaxy is crudely modeled to mimic the pres¬ 
sure confinement by AGN activity, we are confident that such a 
mechanism could operate in more realistic configurations (see the 
jet simulations of Gaibler et al. 2012). We have demonstrated that 
such pressure confinement of the ISM drives the galaxy into an in¬ 
tense star formation regime, and could explain observations of star 
formation-enhanced galaxies in the presence of jet activity (Zinn 
et al. 2013). Cosmological simulations of pure AGN jet feedback 
in galaxy clusters (Dubois et al. 2010) have shown that it has a neg¬ 
ative impact on the galaxy SFR on the long-term, though these sim¬ 
ulations were lacking spatial resolution in order to properly capture 
the small-scale fragmentation of the ISM. Our more global picture 
could suggest a two-stage mechanism for AGN feedback: a com¬ 
pression phase leading to a short burst of star formation, together 
with the expulsion or heating of the circumgalactic gas leading to 
a suppression of the gas accretion onto the galaxy and its star for¬ 
mation on longer time-scales. This remains to be verified with sim¬ 
ulations of galaxies embedded in a cosmological environment with 
high spatial resolution and a self-consistent treatment of AGN feed¬ 
back. We defer this study to future work. 
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APPENDIX A: BIPOLAR PRESSURE INCREASE 

To study the assumption of a isotropic pressure increase, we have 
performed a simulation of a non-isotropic bipolar pressure increase. 
For this the pressure has only been increased after a certain height 
(1.5 kcp) in the vertical direction of the galaxy, were the pressure 
has been kept at the normal value in the radial direction. The SFR 
of the bipolar and isotropic simulations are shown in Fig. Al. One 
can see that while the bipolar SFR oscillates more the general be¬ 
haviour is not changed by the way pressure is applied on the galaxy. 


APPENDIX B: EFFECTS OF SUPERNOVA FEEDBACK 

Here, we compare the feedback run with the no-feedback run. In 
Fig. B1 the gas density maps of the no-pressure enhancement simu¬ 
lations are shown for the non-feedback (nf, left panel) and feedback 
(fb, right panel) simulations. In Fig. B2 the comparison between fb 
and nf is shown for the pa3 simulations. We see that for the no¬ 
pressure simulations, the effect of SN explosions is to disrupt the 
interstellar medium into smaller but more numerous clumps. In the 
edge-on-view, we can also see that the feedback simulation thick¬ 
ens the disc and enhances the mass outflow close to the galaxy. For 
the pressure simulation, no significant difference can be observed. 
It shows that the effect of external pressure is stronger than the ef¬ 
fect of SN explosions. 


gasHigh 



Figure Al. Star formation rate (SFR) as a function of time: In blue for the 
case where the pressure is applied isotropically ( isotropic ) and in red when 
the pressure is applied to the galaxy in a bipolar geometry, (bipolar). 

In Fig. B3, we show the number of clumps as a function of 
time for a selection of the gas Low (left), gasHigh (middle), and 
gasHigh-d (right) simulations with (fb) or without (nf) SN feed¬ 
back. In Fig. B4, we show the time evolution of the SFR for 
the same selection of runs. We see that the number of clumps 
is enhanced by the presence of SN explosions in all cases since 
the clumps are regularly destroyed by the SN activity (Dubois 
et al. 2015). SNe regulate the mass growth of the gas clumps, and 
since the most massive clumps are expected to capture the smaller 
clumps, SNe allow for the increase in the number of clumps, 
thereby reducing their average cross section and mass (see Fig. 7). 
We also see that the SFR is higher for the non-feedback simula¬ 
tion compared to the feedback simulations as a consequence of the 
absence of a local regulating process within gas clumps. 

Reassuringly, the effect of over-pressurisation of the disc onto 
the SFR enhancement is independent of the presence of SN ex¬ 
plosions: it still leads to a positive feedback effect that SNe only 
marginally modulate. 


APPENDIX C: CONVERGENCE STUDIES 

In this section, we test how the results depend on the resolution of 
the simulation. We performed two high resolution ( highRes ) simu¬ 
lations for the gasHigh case, one with no external pressure (nP_hR) 
and the other with external pressure (pa3_hR). The higher resolu¬ 
tion runs have been performed with a resolution of Ax = 10 kpc 
(compared to 40 kpc for the standard runs). We changed the density 
threshold for star formation (no = 224 H cm -3 ) in the poly tropic 
EoS as well as the dissipation time-scale of the non-thermal com¬ 
ponent for the SN feedback (Ax = 10 pc) with the resolution. The 
simulations were run for a similar timescale (~ 400 Myr) as the 
lower resolution ( low Res ) simulations. 

In Fig. Cl, we show the comparison between the highRes and 
lowRes simulations. In the upper panel, the number of clumps is 
shown for the high and low resolution runs where for both simu¬ 
lations the same clump detection density threshold of 21 Hem -3 
and a peak-to-saddle threshold of 1.5 was chosen. 

Fig. Cl shows that, in both highRes and lowRes runs, clumps 
are formed at a faster rate when over-pressure is applied on the 



galaxy. Comparing the two resolution runs, we see that the rates of 
clump formation for both resolutions are comparable at the start of 
the simulations, for both the pressure and no-pressure runs. While 
the lowRes run with external pressure (pa3) sees a sharp rise in its 
clump number at 25 Myr, the number of clumps in the highRes 
run with external pressure (pa3_hR) starts catching up after 50 Myr 
and soon (at 70 Myr) overtakes that of the pa3 run, to end up with 
nearly double the number of clumps. A similar effect is seen in the 
no-pressure runs: the number of clumps in the highRes simulation 
starts slowly, but overtakes that of the lowRes run (at 230 Myr) to 
also end up with nearly double the number of clumps. 

Similar trends are seen in the star formation histories (lower 
panel of Fig. Cl). For the no-pressure runs, the highRes one over¬ 
takes the other one in SFR at 280 Myr to end up with twice the 
SFR, while in the corresponding runs with external pressure, the 
highRes one has its SFR overtake that of the lowRes analog at 150 
Myr, end the highRes run ends up with over double the SFR of the 


lowRes one. The very slow rise of the SFRs in the highRes runs 
is the consequence of our choice of a higher density threshold for 
the highRes simulations, which is reached at later times. Once stars 
start to form, the SFR is greater in the pressure simulation than in 
the no-pressure simulation. The general effect that the pressurisa- 
tion leads to more star formation is therefore still the same. 
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Figure Bl. Gas density maps (mass-weighted) of the gasHigh non¬ 
feedback (left) and feedback (right) simulations without enhancement of 
the external pressure (nP). The maps are taken at the end of the simulation 
(~ 400 Myr). Each panel shows both face-on (40x40 kpc, upper part) and 
edge-on (40x20 kpc, lower part) views. 



Figure B2. Similar as Fig. Bl but for the simulations with pressure en¬ 
hancement pa3. 
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Figure B3. Time evolution of the number of clumps for a selection of the gasLow (left), gasHigh (middle), and gasHigh A (right) gasHigh simulations. For 
each simulation set the feedback (fb) and non-feedback (nf) runs are shown for comparison. They are indicated by the suffixes in the legend. The lines are 
smoothed as in Figure 4. 
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Figure B4. Time evolution of SFR for a selection of the gasLow (left), gasHigh (middle), and gasHighA (right) simulations. For each simulation set the 
feedback (fb) and non-feedback (nf) runs are shown for comparison. They are indicated by suffixes in the legend. The lines are smoothed as in Figure 4. 
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Figure Cl. Time evolution of the number of clumps (upper panel) and 
SFR (lower panel) for the low resolution and high resolution gasHigh-fb 
simulations. The clumps were extracted using a gas density threshold of 
21 H cm -3 . 







